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Two magnetic impurities on the edge of a zigzag graphene nanoribbon strongly interact with each 
other via indirect coupling, which can be mediated by conducting carriers. By means of Quantum 
Monte Carlo (QMC) simulations, we find that the spin-orbit coupling A and the chemical potential 
p in system can be used to drive the transition of local-spin exchange from ferromagnetism to anti¬ 
ferromagnetism. Since the tunable ranges for A and fi in graphene are experimentally reachable, we 
thus open the possibilities for its device application. The symmetry in spatial distribution is broken 
by the vertical and the transversal spin-spin correlations due to the effect of spin-orbit coupling, 
leading to the spatial anisotropy of spin exchange, which distinguish our findings from the case in 
normal Fermi liquid. 
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Due to the unusual electronic quasi-particles and the flexible carrier concentration [lj, graphene with local spins is 
regarded as a novel type of two-dimensional magnetic semiconductor as well as a promising candidate for spintronics 
devices @]. Doping with magnetic impurities is regarded as an effective route to induce local spins in graphene Q, 
because the pristine graphene is a nonmagnetic semimetal without band gap. In previous works concerning single¬ 
impurity doped graphene, it was reported that applied gate voltage (to change the chemical potential) can control the 


electron filling on impurity orbit [4j, [5| and its Kondo screening from conducting electrons [6|, [7[ , thus leading to an on- 
off switching for a single spin in graphene. But in realistic spin-based devices, we need to consider multi-spin case and 
spin-flip process owing to the experimental challenge of doping concentration control, which will induce spin current 
transport in graphene systems. Therefore, for the device application, it is crucial to achieve efficient controllability for 
spin correlations among embedded impurities, which is measured by the Ruderman-Kittel-Kasuya-Yosida (RKKY) 
interaction and mediated by conducting electrons. In previous studies on RKKY interaction, the electron doping 
dominates the large-distance oscillation with a wave vector 2fcp as well as decay behavior I? -7 [§]. 

Recently, there is a surge of interest in effect of spin-orbit coupling (SOC) on electronic properties, which opens 
a class of novel time-reversal invariant topological insulators (TI) [9j. On the edge (two-dimension case) or surface 
(three-dimension case) of these TI materials, the carries have spin-momentum locking, so called ’’helical” states: 
electrons with opposite momenta have orthogonal spin orientations, so that the backscattering k -o- — k is forbidden 
without a spin flip. Due to this scenario, there is spin-current vortex around the single local spin on the helical 
edgeflpj. Furthermore, on the surface of three-dimensional realistic TI, Kondo impurities will open up a local gap 
11] and the exchanges can be described by the RKKY interaction among the impurities [ 13 ]. More interestingly, the 


combination of anisotropic spin exchange and SOC on the edge of two-dimensional realistic TI results in unnormal 
temperature-dependent conductance on the edge [13j, [lj]. These investigations indicate that the SOC can greatly 
modify the spin configuration for single impurity in realistic TI, such as class of Bi 2 Se 3 (Bi 2 Te 3 , Sb 2 Te 3 ) [15| and 
strained HgTe jl6j|. While in graphene, which is one of the earliest proposals for TI materials, multi-spin problem 
has not been comprehensively studied, despite the facts that the SOC in graphene can be enhanced in a wide range, 
while both the Kondo screening around single spin and RKKY interaction between two spins can be greatly affected 
by the SOC. 

In this paper, we consider two Anderson impurities on the edge of a zigzag graphene ribbon, and investigated 
the effect of Kane-Mele SOC and chemical potential on the local spin exchange, which can be achieved by heavy 
metal doping or electric field. Through Quantum Monte Carlo simulation at finite temperature, we carry out a 
comprehensive study on the spin and charge coupling strengthes of this impurity pair. For the spin-spin correlation, 
which measures the RKKY interaction, we found that, the in-plane components of correlations favor ferromagnetism 
but the out-plane correlation can be tuned from ferromagnetism to antiferromagnetism by the SOC, while such 
transition is much dependent on the on-site Coulomb interaction of impurity orbital. From recent first-principal and 
experimental results of enhanced SOC in graphene, we conclude that spin-spin exchange among impurities in graphene 
is experimentally tunable by SOC. In contrast, the charge fluctuation between two impurities is always suppressed. 
Meanwhile, in-plane components of correlation will decay with an oscillation behavior when the distance between two 
impurities is increased, but the out-plane spin and charge correlations will decay smoothly, which is totally different 
from the Fermi liquid behavior, and leading to the spatial anisotropy of spin exchange. 
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To simulate the model, we construct a model with width of 20 A zigzag-edge graphene ribbon and along the 
horizontal direction we adopt periodical condition. As indicated in Kane-Melc investigation [17], graphene will be 
driven into topological phase when the intrinsic SOC in graphene is enhanced. In the case, the ribbon would possess 
the helical edge states. In Fig. [I] we show the scheme for two impurities (yellow balls) located on the zigzag helical 
edge, where red and blue arrows denote the conducting electrons with up and down spins, respectively, and they 
propagate along oppositive directions (spin-momentum locking). 

The QMC simulation is based on two-Anderson-impurity model, the total Hamiltonian reads as 


H = H k _ m +H 1 +H 2 , (1) 

where F/k-m is the Kane-Mele model in a zigzag edge graphene ribbon il7|. It is composed with two parts: Hk~m = 
H t + H so , with H t being the usual nearest-neighbor hopping of tight-binding model in graphene 

Ht = —t ^ ' ( C i,cr C i,o’ T H.C.) — M ^ ^ c i.a c i,&’ 

<ij>,a i,a 


and 


H so = X ^2 (ivijc\(j z Cj + H.c.) 

«ij» 

H so is the SOC term and a z is the z Pauli matrix. H so thus has opposite signs for opposite electron spins. The 
parameters v,jj = —Vji = ±1 depend on the orientation of the two nearest neighbor bonds as the electron hops from 
site i to j: Vij = +1 if the electron makes a left turn to the second bond. It is negative if it makes a right turn. H\ 
is the part for two impurities 


— 'y ' ( £ ™ AO^crm^o-m + U m d^ m d^ m d^ m d^ m . 


<7771=1,2 

and if 2 describes the hybridization between impurities and conducting electrons 


FF 2 = 


\ ' 


2 — 

am— 1,2 


V m [ C ma^crm A Arm Oner] . 


Here m is the index of two impurities, in Hi and 7F 2 , e m is the energy of the impurity orbital and U m is Coulomb 
repulsion inhibiting the simultaneous occupancy of the orbital by two electrons and V m is the hybridization strength. 
In this work, we consider the two identical impurities, so the above parameters in Hamiltonian are m-independent. 

We use the single-impurity QMC base on discrete-time algorithm [18J to compute the local thermodynamic prop¬ 
erties of the impurity. The length of imaginary-time is set as At = 0.1 to suppress the systemic error. The QMC 
naturally returns the imaginary-time Green’s function Gd<j{r > 0) = (d<j(r)dJ.) of the impurity. With this Green’s 
function, we can compute the spin-spin correlations for two impurities along a = x,y,z directions versus distance as 


S a (R) = <«>, 

where s x m = ( d\ m d im + d\ m d fm ), s v m = {-i) 2 (d\ m d lm - d\ m d tm ) and s z m = {d\ m d tm - d\ m d im )-, as well as the 
charge-charge correlation as the function of distance 


C{R) = {nin 2 ) - (ni)(n 2 ), 

and n m = (dld^rn + d\ m d, m ). The spin-spin correlations S a (R) directly reflect the RKKY interaction which is 
mediated by conducting electrons on the helical edge, while charge-charge correlation C(R) indicates the indirect 
Coulomb repulsion. 

In Fig. [2j we show the spin-spin and charge-charge correlations as function of SOC, Coulomb interaction, and 
chemical potential. Here, we fixed the hole-particle symmetry n = 0 and e = — y, so that {di m d^ m + d^ m d^ m ) = 1 
and the local moments at two impurity sites are well developed. One can see from Fig.[2](a), the out-plane correlations 
of S z highly depend on SOC, which can be reduced in a linearly way and have a wide variation range. When SOC is 
increased, the value of vertical correlation S z can be changed from positive to negative, it means that a transition from 
a ferromagnetic (FM) to anti-ferromagnetic (AFM) correlation is achieved when the SOC reaches a critical value. It is 
interesting to note that, the critical SOC strength for such a transition increases as Coulomb interaction is increased, 
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from the QMC simulations. For example, the critical SOC value is about 100 meV when U is leV, but it is increased 
to 300 meV when U is is 5eV. The underlying mechanism can be attributed to the competition of spin and charge 
fluctuations between two impurities. When they are close to each other, the changing in carriers SOC can directly 
result in the spin-flip process which reflects in S' 2 , but the increasing Coulomb interaction U can block spin flip, so 
it is seen that with a larger U in Fig. HJa), the conducting electrons need larger kinetic energies (SOC) to drive a 
FM-AFM transition. But when two impurities are far away from with each other, we can not see this transition in 
S z but in S x ( y \ which is totally determined by the Luttinger liquid physics on the edge [Tol ]. The explanation can be 
also demonstrated in the latter plots [Fig. 3(a)]. 

The SOC driven magnetism transition will be important for the spintronics application. It is worthy to point out 
the required critical SOC is in a reasonable range based on the fact that recent extensive strategies were proposed 
and developed to enhance SOC in graphene. In Table [I] we sumarize the recent achievements from both density 
functional theory and experiments. In the case of interface/intercalation, the recent tunable range for SOC is 50 - 
225 meV through placing graphene on the substrates of Au(lll), Ag(lll), and Ni(lll) and intercalation with SiC 
(oooi)Eaa Alternatively, doping graphene with heavy or magnetic adatoms (Ni, Au, Ti and Os) is also able to 
induced the enhanced SOC in graphene of 20 - 500 meV [25l427| . Meanwhile, when sandwiching graphene between 
strong spin orbit interactions or Rashba materials by the proximity effects, the intrinsic SOC is significantly enhanced 
and graphene is driven into topological phase with nontrivial gap of 120 meV mm. Due to these efficient strategies, 
it is possible to tune the SOC in graphene over a wide range, the FM-AFM transition for spin correlation between 
two magnetic impurities via SOC controlling is thus feasible. 

In contrast to the spin-spin correction transition of S z , the transversal correlation S x ^ y ' ) increased linearly by the 
SOC interaction but not change its sign. In Fig. [2](b), we see that at SOC A = 0, S z = S x ' v \ but this symmetry is 
broken when SOC is switched on. The transversal correlation S X<JJ> is increased when either of SOC and Coulomb 
interaction is enhanced. More interestingly, the enhance SOC not only can change the spin correction between two 
impurities, but also affect the charge interaction between them. As indicated in Fig.[2](c), the charge-charge correlation 
as a function of SOC with different values of Coulomb interaction is calculated. We can see that the absolute value of 
C is increased when SOC is increased, but decreased as the Coulomb interaction is enhanced. When U is at 1.2 1, the 
charge fluctuation is totally suppressed. To have a comprehensive physics picture, we also vary chemical potential ^ 
to break the hole-particle symmetry and study their effects on S z and S x ^ v \ In Fig.[2jd). The vertical and transversal 
correlations are totally degenerated (i.e. S z is equal to S x ^) at SOC strength A = 0. But when A ^ 0, the SOC effect 
on both spin correlation is different, leading to non-degenerated S z and S x ^ v \ We then turn to the effect of chemical 
potential on S z and S x ^ v \ it could be found that both S z and S x< ' v ' ) are positive near the zero point of /it, suggesting 
that two impurities correlate ferromagnetically. After /i is shifted from 0, the two impurities have anti-ferromagnetic 
correlation as indicated by the negative values of S z and S x l y ) ( ^ is in the range of -0.2t to -0.8t). But when is 
further reduced, both S z and approach to zero, indicating the weak coupling between two impurities. 

The distance between two impurities is expected to influence the couplings of spin-spin and charge-charge due to 
the fact that the RKKY interation is sensitive to the distance. We preformed the QMC calculations to address the 
problem when varying the distance. In Fig. [3] we show the spin-spin and charge-charge correlations as functions of 
distance between two impurities. Here we consider the case that both of two impurities are located on the edge (the 
same sub-lattice). In Fig. (3][a), we see that the vertical spin-spin correlation S Z (R) and transversal one S x ( y \R) 
have different spatial distribution. S x( - V \R) has both the behavior of decaying and oscillating, but S Z (R) is always 
decaying with the increased distance. Similar to S Z (R), the charge-charge correlation is also decaying with distance 
(Fig. EKb)). These correlation functions have the following behavior: 


S Z (R) oc 


h(RY 


S-M(R) oc 

h{R) 


and 


C{R) oc 


fs(R)’ 


here fi(R), f 2 {R) and fa(R) are decaying function of distance has the formula as l/R?. 

To more clearly display these correlations, we calculate the ratios S X ^ V \R)/C(R) and S Z (R)/C(R). We find that 
S Z (R)/C(R) is a trivial constant, but S x ^ y \R)/C(R) is an oscillating function of distance, as shown in Fig.[3](c) and 
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(d). The ratio exhibits prefect behavior of sine or cosine functions with wave vector, which depends on the chemical 
potential and SOC. 

Finally, we discuss the model Hamiltonian of Eq. m in more realistic situations with the effect of Hubbard U' in 
graphene and hybridization. In realistic graphene zigzag nanoribbon, the Coulomb interaction plays a signification 
role and stabilizes a spin-polarized ground state [301. In this case, the conducting electron at the edge can be regarded 
as standard one-dimensional Luttinger liquid, coupled with which the two magnetic impurities at low temperature 
have the RKKY interaction as the asymptotic behavior I ? -7 cos 2 kpR, here kp is the Fermi momentum. The decaying 
factor 7 is the function of Hubbard IP: when U' = 0, 7 = 1, and when U' ^ 0, 0 < 7 < 1 [31, HJ. So in Luttinger 


liquid case, at low temperatures the RKKY interaction decays more slowly than the case in Fermi liquid U' = 0 . 
The methods used to enhance the graphene’s SOC in Table Q] are based in the mechanism that via the hybridization 
between the electronic bands of carbon atoms and heavy atoms with large intrinsic SOC, the SOC of graphene will be 
greatly enhanced. The hybridization depends on three parameters of heavy atoms: the crystal field, the intrinsic SOC 
and the hopping between carbon and external atoms. If the crystal field is greatly larger than the later two quantities, 
the graphene 7 r bands still dominate near the Fermi energy and Kane-Mele Hamiltonian is also valid and the effect 
of hybridization can be absorbed in the enhanced SOC through second order approximation [2(| and furthermore 
when the crystal field decreases and the intrinsic SOC and the hopping between carbon and external atoms increase, 
the SOC in graphene will increase. While if these three parameters have the same order and crystal field is close 
to Dirac point, on one hand such strong hybridization can further enhance the SOC in graphene and influence the 
RKKY interaction as shown in Fig. [2ja)-(c); on the other hand, there are additional bands near the Fermi energy 
from heavy atoms, the electrons from these bands will exchange with those at the levels of impurities (both spin and 
charge exchanging). As a result of this exchanging, the spin-spin and charge-charge correlations will increase. In 
addition, the three ways in Table U can raise charge transferring and shift the Fermi energy of pristine graphene from 
Dirac point, so we do the calculation varying the chemical potential and show the effects in Fig. [2])d) . 

In conclusion, we investigate two Anderson impurities on a zigzag graphene nanoribbon by QMC simulation in 
present work. When the system is under a topological phase transition driven by Kane-Mele SOC, we find that , the 
spin rotation symmetry in RKKY interaction between two impurities can be partly broken because of the SOC. As the 
two impurities are close to each other and the Fermi energy is near the Dirac point, the in-plane RKKY interactions 
favor ferromagnetism but the out-plane component can be tuned from ferromagnetism to anti-ferromagnetism by the 
SOC. The distance effect on RKKY interaction is also systematically investigated, the in-plane interaction has both 
behavior of decaying and oscillation, but the out-plane interaction and charge-charge correlation only has decaying 
as the distance between them is increased. This behavior is totally different from the conventional Fermi liquid. 
Our theoretical findings here provide some useful suggestions of the spin-polarized STM measurement to RKKY 
interaction|33i|. 
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FIG. 1: (color online) The scheme for the two magnetic impurities (yellow balls) at the edge of the zigzag nanoribbon. The 
arrows on the top and bottom demotes the moving directions for two species spins. 



FIG. 2: (color online) (a) The spin-spin correlation along z direction S z versus A, (b) The spin-spin correlation along x or 
y direction S x ' y ^ versus A (c) The charge-charge correlation versus A and in (a)-(c) the inverse of temperature /3 = 321 _1 , 
V = 0.65 1, fi = 0 and e = — y. (d) S z and S x ^ versus chemical potential n with /3 = 32t _1 , U = 1.2f, V = 0.65 1 and e = —y. 



FIG. 3: (color online) (a) The vertical and transversal spin-spin correlation as function of distance between two impurities (b) 
The charge-charge correlation as function of distance, (c) The ratio S X< ' V \R)/C(R) with fixed chemical potential /r = — 0.2t 
and varied SOC A = O.lf and 0.15t. (d) The ratio S X< ' V \R)/C(R) with fixed SOC A = O.lt and varied chemical potential 
fj, = O.lf and 0.2f.The x axis is in the unit of lattice constant of graphene. 





















































7 


TABLE I: Enhanced Gap 

AEsoc 

in Graphene 

Method 

Asoc (meV) 

Interface/intercalation 

50 - 

225 

[20-24] 

Adatom doping 

20 - 

500 

[25-27] 

Proximity effect 

50 - 

120 

[28, 29] 








